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ABSTRACT 

Gaudi & Gould (1997) showed that close companions of remote binary systems can 
be efficiently detected by using gravitational microlensing via the deviations in the 
lensing light curves induced by the existence of the lens companions. In this paper, 
we introduce another channel to detect faint close-in binary companions by using 
microlensing. This method utilizes a caustic-crossing binary lens event with a source 
also composed of binary stars, where the companion is a faint star. Detection of the 
companion is possible because the flux of the companion can be highly amplified when 
it crosses the lens caustic. The detection is facilitated since the companion is more 
amplified than the primary because it, in general, has a smaller size than the primary, 
and thus experiences less finite source effect. The method is extension of the previous 
one suggested to detect close-in giant planets by Graff & Gaudi (2000) and Lewis & 
Ibata (2000) and further developed by Ashton & Lewis (2001). From the simulations 
of realistic Galactic bulge events, we find that companions of K-type main sequence 
or brighter can be efficiently detected from the current type microlensing followup 
observations by using the proposed method. We also find that compared to the method 
of detecting lens companions for which the efficiency drops significantly for binaries 
with separations < 0.2 of the angular Einstein ring radius, #e, the proposed method 
has an important advantage of being able to detect companions with substantially 
smaller separations down to ~ (0)1O _2 0e- 
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1 INTRODUCTION 

Observationally constraining the characteristics of faint 
close-in binary companions is important to understand the 
process of star formation and evolution. For example, mass 
transfer between the binary components through Roche lobe 
overflow could alter the physical processes of stellar evolu- 
tion (Iben & Tutukov 1984; Han, Tout, & Eggleton 2000). 
Theoretical investigations about the origin of Type la su- 
pernova are also based on the observations of close bina- 
ries (Nomoto & Sugimoto 1977; Li & van den Heuvel 1997; 
Langer et al. 2000). Another possible outcome of close bi- 
nary evolution is a merger during the late stage of a common 
envelope phase (Iben & Livio 1993; Armitage & Livio 2000; 
Taam & Sandquist 2000) , whose physical parameters should 
be constrained by observations. 



Detection of faint close-in binary companions is also im- 
portant for the accurate determination of the stellar mass 
function. The mass function of stars is derived from their lu- 
minosity function by using the mass-luminosity relation. In 
this determination, however, if the mass function is derived 
without considering unresolved faint binary companions, the 
resulting stellar mass function will be biased towards larger 
masses. Even if they are taken into consideration, the con- 
structed mass function will suffer from large uncertainties 
because the uncertainties in the properties (the frequency 
and the mass ratio distribution) of the unresolved bina- 
ries will directly propagate into the resulting mass function 
(Kroupa, Tout, & Gilmore 1991; Reid 1991; Reid & Gizis 
1997). 

Faint close-in binary companions can be detected by 
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using various methods, e.g. direct imaging by using high- 
resolution interferometers, astrometric measurements of the 
primary motion, spectroscopic radial velocity measure- 
ments, and photometric measurements of the light varia- 
tion caused by occultation (eclipsing) . However, these meth- 
ods have limitations in detecting and characterizing various 
types of binaries. For example, the direct imaging and astro- 
metric methods have limited applicabilities only to nearby 
stars. The spectroscopic method can be applied to detect 
binaries located at large distances. However, determining 
the mass ratio between the binary components by using 
this method requires that both components should be bright 
for the construction of the individual stars' radial velocity 
curves, limiting the applicability of the method. Constrain- 
ing the physical properties of eclipsing binary companions 
suffers from similar limitation. Therefore, invention of com- 
plementary methods to detect and characterize various types 
of binaries is important. 

Close binaries can also be detected and characterized 
by using microlensing. Detecting binary companions via mi- 
crolensing is possible because if events are caused by lenses 
composed of binaries, the resulting light curves can exhibit 
noticeable deviations from the smooth and symmetric ones 
of single lens events. In addition, one can constrain the phys- 
ical properties of the binaries (the mass ratios and the scaled 
projected separation between the binary components) by fit- 
ting model light curves to the observed one. 

In this paper, we introduce another method of detect- 
ing and characterizing faint close-in binary companions by 
using microlensing. This method utilizes caustic-crossing bi- 
nary lens events with a source also composed of binary stars, 
where the companion is a faint star. Detecting the source 
companion is possible because the flux from the companion 
can be highly amplified when it crosses the lens caustics. 
This method is extension of the previous one suggested to 
detect close-in giant planets by Graff & Gaudi (2000) and 
Lewis & Ibata (2000) and further developed by Ashton & 
Lewis (2001). 

The paper is organized in the following way. We de- 
scribe the basics of microlensing in § 2. We perform simula- 
tions of example Galactic bulge events in § 3 and evaluate 
the detectability of the deviations induced by binary source 
companions for the example events in § 4. A brief summary 
of the results and discussion about the advantages and dis- 
advantages of the proposed method compared to that of 
detecting lens companions are included in § 5. 



2 BASICS OF MICROLENSING 

The light curve of a single point-source event caused by a 
single point-mass lens is represented by a simple analytic 
form of 



A = 



u^u 2 + 4 



(1) 



where u is the lens-source separation normalized by the an- 
gular Einstein ring radius. The Einstein ring radius is related 
to the physical parameters of the lens by 
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(2) 



where M is the mass of the lens and D a \ and D os repre- 
sent the distances from the observer to the lens and source, 
respectively. 

The light curve of a caustic-crossing binary lens event 
with a source composed of also binaries, which is the target 
event for followup monitoring to detect faint binary source 
companions by using the proposed method, deviates from 
the standard form in eq. (1) due to various effects. These 
include the effects of lens and source binarity, and the finite 
source effect. In the following subsections, we describe the 
basics of these effects. 



2.1 Binary Lens Events 

If an event is caused by a binary lens system, the locations 
of the individual images are obtained by solving the lens 
equation expressed in complex notations by 



nil 



+ 



III > 



(3) 

Z\ — Z Z2 — Z 

where mi and rri2 are the mass fractions of the individual 
lenses, z\ and Z2 are the positions of the lenses, f = £ + ir] 
and z = x + iy are the positions of the source and images, 
and z denotes the complex conjugate of z (Witt 1990). Note 
that all these lengths are dimensionless because they are 
normalized by the combined Einstein ring radius, which is 
equivalent to the Einstein ring radius of a single lens with a 
mass equal to the total mass of the binary. Depending on the 
source position with respect to the lenses, there exist 3 or 5 
solutions for the lens equation, and thus the same number 
of images. The amplification of each image A; is given by 
the Jacobian of the lens equation evaluated at the image 
position, i.e., 



A ,i = 



|det J\ 



az az 



(4) 



Then, the total amplification is obtained by summing the 
amplifications of the individual images, i.e. Aq = ^ Ao,-;. 

The fundamental difference in the geometry of a binary 
lens system from that of a single point-mass lens is the for- 
mation of caustics. The caustic refers to the source position 
on which the amplification of a point source event becomes 
infinity, i.e., det J — 0. The set of caustics forms closed 
curves. The caustics takes various sizes and shapes depend- 
ing on the projected separation b (normalized by #e) and 
the mass ratio q between the lens components. The proba- 
bility of caustic crossing is maximized when the binary sep- 
aration is equivalent to 9e- Whenever a source crosses the 
caustics, an extra pair of images appears (or disappears). 
Hence the light curve of a caustic-crossing binary lens event 
is characterized by sharp spikes. Since the caustics form a 
closed curve, the source of a caustic-crossing event crosses 
the caustic at least twice. Although the first caustic crossing 
is unlikely to be resolved due to its unpredictable and swift 
passage, it can be inferred from the enhanced amplification. 
Then, by preparing intensive followup observations, resolv- 
ing the second caustic crossing is possible. From the surveys 
(EROS: Aubourg et al. 1993; MACHO: Alcock et al. 1993; 
OGLE: Udalski et al. 1993; DUO: Alard & Guibert 1997) 
and followup observations (MPS: Rhie et al. 1999; PLANET: 
Albrow et al. 1998; MOA: Abe et al. 1997), numerous binary 
lens event candidates have been detected and a significant 
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fraction of them are caustic-crossing events (Alcock et al. 
2000). 

2.2 Finite Source Effect 

In most parts, the light curve of an observed caustic-crossing 
binary lens event is well described by the point-source ap- 
proximation. However, the part of the light curve near and 
during the caustic crossing deviates from this approxima- 
tion. This is because the source star in reality is not a point 
source and the gradient of the amplification near the caustics 
is very large. As a result, different parts of the source star 
disk are amplified by considerably different amounts even 
with small differences in distance to the caustics (Schneider 
& Weiss 1986; Witt & Mao 1994). 

The light curve of an extended source event is given 
by the intensity-weighted amplification averaged over the 
source star disk, i.e., 

= JTX I(r,e)M\r-r L \)rdr d9 
frS? I{r,6)r dr dO 

where r* is the radius of the source star, 7(r, 9) is the surface 
intensity distribution of the source star, and the vectors r 
and vl represent the displacement vector of the center of 
the source star with respect to the lens and the orientation 
vector of a point on the source star surface with respect 
to the center of the source star, respectively. Consequently, 
the amplification does not become infinite even during the 
caustic crossing. 

2.3 Binary Source Events 

The light curve of an event with a source composed of double 
stars (binary source event) is very simple, as it is just the 
linear sum of the two light curves of the individual single 
source events (Griest & Hu 1992). Since the baseline flux is 
also the sum of the (unamplified) fluxes from the two source 
stars, the amplification becomes 

F ,i + Fo, 2 v ; 

where Fo,i and Ai are the baseline fluxes and the amplifica- 
tions of the individual single source events and the subscripts 
i=l and 2 denote the primary and the companion source 
stars, respectively. 



3 SIMULATIONS 

The detectability of the light curve deviations induced by 
close-in binary source companions depends on various com- 
binations of many parameters such as those characterizing 
the binary source system (e.g., the projected separation and 
the relative phase between the source components) and the 
individual component source stars (e.g., the brightness and 
size). Consequently, it is difficult to present the detection 
efficiency as a function of a couple of parameters. There- 
fore, we choose to evaluate the detectability by testing ex- 
ample light curves produced from the simulations of typical 
Galactic bulge events. The simulations are performed in the 
following ways. 



Table 1. The adopted values of the physical parameters of the 
source stars for the tested example lensing events. The values 
are adopted from Allen (1973). 



stellar 


absolute 


radius 


type 


magnitude 


(«©) 


Clump Giant 


+2.9 


4.0 


GO Dwarf 


+4.4 


1.04 


K0 Dwarf 


+5.9 


0.85 


MO Dwarf 


+9.0 


0.63 



For the source stars, we test various combinations of 
stellar types, which well characterize the individual source 
components. In Table 1, we list the source stellar types of the 
tested events along with the absolute magnitudes and radii, 
which are adopted from Allen (1973). In the simulation, we 
assume that each source star has a uniform surface bright- 
ness distribution over its diskQ and it is located at D os = 8 
kpc. We also assume that the motion of each source is rec- 
tilinear. This is not only because the changes of the source 
positions caused by their orbital motion will be small during 
the caustic crossings, but also because what directly affects 
the detectability are the instantaneous projected separation, 
I (normalized by 9e), and the phase angle, a, between the 
source components at the time of caustic crossing. 

For the lens system, we assume that the mass ratio and 
the projected separation between the lens components are 
q — 0.5 and b = 1.2. The assumed Einstein ring radius 
crossing timescale (Einstein timescale) of the events is 4b = 
27 days, which corresponds to the value of a Galactic bulge 
event caused by a lens with a total mass and location of 
~ 0.3 Mq and D \ ~ 6 kpc under the assumed lens-source 
transverse speed of v ~ 150 km s _1 . 

In Figure 1, we present the caustics of the lens sys- 
tem and the source trajectories of one of the tested example 
events. We note that other tested events have similar tra- 
jectories except the variations caused by the differences in 
the values of £ and a. 

fn Figures 2 and 3, we present the light curves of events 
resulting from the various combinations of source stellar 
types. We present only the part of the light curve around the 
second caustic crossing of the primary, because high time- 
resolution observation is possible at around this time. For 
the light curves in Fig. 2, the primary source star is a red 
clump giant (RG) and the corresponding companions are 
GO, K0, and M0 main-sequence (MS) stars, respectively. For 
the light curves in Fig. 3, the primary is a GO main-sequence 
and the companions are K0 and M0 main-sequence stars, re- 
spectively. Note that the two tested stellar types of the pri- 
mary source are the typical ones of Galactic bulge events. 
For these events, we assume a common separation and a 
phase angle of I — 0.1 and a = 160°. Since the Einstein ring 



* If the source star disk is not uniform due to, for example, limb 
darkening (Witt 1995; Loeb & Sasselov 1995; Gould & Welch 
1997; Valls-Gabaud 1998; Gaudi & Gould 1999 ; Ignace & Hendry 
1999) or spots (Hcyrovsky & Sasselov 2000; Han et al. 2000), the 
resulting light curve deviates from the uniform disk approxima- 
tion. However, we note that the amount of this additional devia- 
tion is slight. 
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Figure 1. The location of the caustics (the closed figure drawn 
by a thick solid line) and the source star trajectories (straight 
lines) of one of the tested caustic-crossing binary lens events with 
double source stars. The coordinates are chosen so that the center 
of mass of the binary lens is at the origin. Both lenses are on the 
£ axis and the heavier lens is to the right. The separation and the 
mass ratio between the binary lens components are b = 1.2 and 
q = 0.5. The inset shows the locations of the binary source stars 
and the phase angle at a particular moment. 



radius projected on the source plane of the example events 
is f-E = r-ED os / D \ ~ 2.6 AU, the assumed source separa- 
tion corresponds to ~ 0.25 AU or ~ 56 Rq. We define a 
as the angle between the primary source trajectory and the 
line connecting the primary and secondary source stars mea- 
sured in the clockwise sense (see the inset in Fig. 1). From 
the figures, one finds that the part of the light curve dur- 
ing the caustic crossing of the secondary source is spikier 
than the part during the primary source caustic crossing, 
implying that the secondary is more highly amplified. This 
is because the secondary, in general, has a smaller size than 
the primary, and thus experiences less finite source effect. 

To investigate the dependency of the detectability on 
the binary source separation, we also perform simulations 
of events occurred on a binary source stars with various 
values of £, and the resulting light curves are presented in 
Figure 4. For these events, the stellar types of the source 
stars are a red clump giant for the primary and a K0 main 
sequence for the companion, respectively. We test four cases 
of source separation of £ =0.01, 0.05, 0.1, and 0.15. To see 



a 3 



2 - 



1 - 



RG (single source) 
RG + GO MS 
RG + K0 MS 
RG + M0 MS 




-0.05 



0.05 



Figure 2. The light curves of a caustic-crossing binary lens event 
with a source also composed of binary stars where the primary is 
a red clump giant. The curves with different line types represent 
those expected when the primary source is accompanied by com- 
panions of GO, K0, and M0 main-sequence stars, respectively. The 
time is normalized by the Einstein time scale and it is measured 
from the moment when the center of the primary star crosses the 
caustic. The source stars have a common separation and a phase 
angle of I = 0.1 and a = 160°. 



the dependency of the detectability only on the separation, 
we assume a = 180°, implying that both stars are aligned 
and move along a common trajectory. 

We also investigate the dependency of the detectability 
on the relative source phase by producing light curves of 
events occurred on source stars with various phase angles. 
In Figure 5, we present the resulting light curves. We test five 
cases of the phase angle of a = 0°, 45°, 90°, 135°, and 180°. 
For these events, we also assume a red clump giant primary 
and a K-type main sequence companion. The source stars 
have a common separation of I = 0.1. 



4 DETECTABILITY EVALUATION 

In this section, we evaluate the detectability of the binary 
source induced deviations for the events produced by simu- 
lations in the previous section. 

First, to be readily noticed, the light curve should de- 
viate from that of a single source event with a significant 
fractional deviation. The major part of the binary source 
event light curve is well approximated by that of a single 
source event with a baseline source flux equal to the sum of 
the fluxes from both source stars. We, therefore, define the 
fractional excess amplification by 

A -A s 



(7) 



where A and A s represent the amplifications of the binary 
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Figure 3. Similar light curves as those in Fig. 2, but for events 
with different combinations of source stellar types. The individual 
light curves are those expected when the primary of GO main- 
sequence is accompanied by companions of KO and MO main- 
sequence stars, respectively. The source stars have a common 
separation and a phase angle of I = 0.1 and a = 160°. 



5 - 1 = 0.15 (RG + KO MS, a=180°) 

; 1=0.1 ; 

1 = 0.05 




0.05 0.1 0.15 

t 



Figure 4. Light curves of binary lens events occurred on binary 
source stars with various projected separations. The stellar types 
of the source stars are a red clump giant for the primary and a K0 
main sequence for the companion. The assumed source orientation 
angle is a = 180°. 



Figure 5. Light curves of binary lens events occurred on binary 
source stars with various phase angles. The assumed stellar types 
of the primary and companion are a red clump giant and a K- 
type main sequence, respectively. The source stars have a common 
separation of i = 0.1. 

and single source events, respectively. Then, we assume that 
the source companion can be detected if the fractional de- 
viation is greater than a threshold value of e t h- 

Second, even if the fractional deviation is substantial, 
the deviation cannot be detected if the absolute amount of 
the deviation is smaller than the photometric uncertainty. 
Currently, microlensing followup observations are being (or 
planned to be) carried out by using the difference image 
analysis method for better photometric precision (D. Ben- 
nett, private communication). This method measures the 
light variation by subtracting an observed image from a 
convolved and normalized reference image of the same field 
(Tomaney & Crotts 1996; Alard & Lupton 1998). Then, the 
signal measured on the subtracted image is proportional to 
the source flux variation, i.e. S oc Fo(A — 1), while the noise 
originates from both the source and background flux, i.e. 
N (X FqA + B, and thus the resulting signal-to-noise ratio is 
approximated by 

where Fo = i*o,i +fo,2 is the total baseline flux of the source 
stars, tcxp = 2 min is the adopted exposure time, and (B) is 
the average background flux within the point spread func- 
tion (PSF). We assume that the observations are carried 
out in I bandQ by using aim telescope with a CCD camera 
that detects 12 e~ s _1 for a star with / = 20 mag. Terndrup 

t This is because in / band more photons can be detected due to 
its broad band width and small amount of extinction. 
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Table 2. The durations of deviations induced by the binary source compan- 
ions for the tested example lensing events. For the definition of td U r, see § 
4. 







source stars 




^dur 


related 


stellar type 


separation (#e) 


phase angle 


(hrs) 


figure 


RG + 


GO MS 


0.1 


160° 


46.0 


Fig. 2 


RG + 


KO MS 






32.4 




RG + 


MO MS 


_ 


_ 


0.0 


_ 


GO MS 


+ KO MS 






46.0 


Fig. 3 


GO MS 


+ MO MS 






0.0 




RG + 


KO MS 


0.01 


180° 


4.5 


Fig. 4 






0.05 




28.5 








0.1 




28.5 








0.15 




28.5 




RG + 


KO MS 


0.1 


0° 


1.5 


Fig. 5 








45° 


7.8 










90° 


5.8 










135° 


5.2 










180° 


28.5 





(1988) determined that the dereddened background surface 
brightness contributed by bulge stars is hi ~ 17.6 mag 
arcsec -2 . With the adopted mean extinction of {Ai) = 1.5 
towards the Galactic bulge field (Stanek 1996), and ad- 
ditional 30% background flux from disk stars (Terndrup 
1988) and sky background flux contribution of ~ 1/4 of 
stellar background, the average background flux within the 
PSF (with an angular area of fi ps f ~ 3.1 arcsec 2 ) becomes 
(B) ~ 560 e _1 s _1 . Once the signal-to-noise ratio is com- 
puted, the uncertainty of the light variation measurement in 
magnitude is computed by 

_ SF/F 1.09 

0.4 In 10 ~ S/A" 1 ' 

where SF/F is the fractional uncertainty in the measured 
flux. Then, if the deviation occurs when A ~ 2, the signal- 
to-noise ratio becomes S/N ~ 14 for the event with the 
primary source of the clump giant, and thus the photometric 
uncertainty is a ~ 0.08 mag. For the event with the primary 
source of the GO main-sequence, the signal-to-noise ratio is 
S/N ~ 3.7 and the corresponding photometric uncertainty 
is a ~ 0.29 mag. 

The detection is additionally restricted by the duration 
of the deviations, tdur- We, therefore, estimate tdur for the 
example events by defining td ur as the length of time period 
during which the fractional and absolute deviations greater 
than e t h = 0.1 and 3<r level, respectively. The resulting val- 
ues of tdur of the individual example events are listed in 
Table 2. Considering the monitoring frequency of the cur- 
rent followup observations of ~ 2 hr _1 , we assume that the 
deviation can be detected if it continues for tdur > 5 hr. 

The findings from the evaluation of the detectability are 
summarized as follows: 

(i) For both red clump giant and G-type main-sequence 
primaries, we find that deviations induced by source com- 
panions brighter than K-type main-sequence or brighter can 
be efficiently detected by using the proposed method. How- 



ever, detecting deviations induced by M-type main-sequence 
companions will be difficult. 

(ii) In addition, we find that with the proposed method 
it will be possible to detect companions with separations as 
small as I ~ 0.01. For a typical event with f e ~ 2 - 3 AU, 
this separation corresponds to ~ 5 Rq. Considering that the 
size of the a clump giant is ~ 4 Rq, this implies that one 
can detect companions of nearly contact binary systems.^] 
We also find that as long as the separation is larger than 
the limiting value, the dependency of the detectability on 
the source separation is weak. 

(iii) Unlike the weak dependency on I, we find that the de- 
tectability depends strongly on a: it will be easier to detect 
the deviation if it occurs after the primary caustic crossing. 
Although it will be possible to detect the deviation occurred 
before the primary caustic crossing, the duration of the de- 
viation is significantly shorter than that of the deviation 
occurred after the primary caustic crossing. This is because 
the amplification of the primary source flux at the time of 
deviation is high, and thus the fractional deviation caused 
by the companion is low. 



5 SUMMARY AND DISCUSSION 

In this paper, we propose a new method of detecting 
faint close-in stellar companions by using microlensing. The 
method utilizes a caustic-crossing binary lens event with a 
source also composed of binary stars, in which the compan- 
ion is a faint star. The detection is possible because the 
flux from the companion can be highly amplified when it 

t The light curve in such cases will be different from that pre- 
sented in Fig. 4 due to the effect of stellar distortion. However, 
analyzing the light curve variations caused by this effect is be- 
yond the scope of this paper, and thus we use the approximation 
of circular source shape. 



© 0000 RAS, MNRAS 000, |-?? 



Microlensing Faint Binary Companion Detection 7 



crosses the lens caustic and is facilitated because the com- 
panion is likely to be more amplified than the primary due 
to the smaller size of the companion. From the simulations 
of Galactic bulge events under realistic conditions, we find 
that companions with K-type main-sequence or brighter of 
binary lens systems with separations down to ~ (0)1O -2 
of the angular Einstein ring radius can be effectively de- 
tected by using the proposed method from the current type 
followup observations. 

The proposed method of detecting source companions 
has both advantages and disadvantages compared to that 
of detecting lens companions. The most important disad- 
vantage is that very faint or dark companions cannot be 
detected by using the proposed method, while the latter 
method can detect companions regardless of their luminosi- 
ties. However, the proposed method has an important ad- 
vantage of being able to detect very close companions with 
separations as small as ~ (0)1O _2 #e- By contrast, the effi- 
ciency of the latter method drops rapidly for binaries with 
separations < 0.2# B (Gaudi & Gould 1997). Therefore, the 
proposed method will make microlensing an important tool 
to detect binary companions by complementing the previous 
one. In addition, since one can obtain the color and spec- 
tral information of the companion from followup multi-band 
photometric and spectroscopic observations performed dur- 
ing the deviation by taking advantage of its highly amplified 
flux, one can characterize the stellar type of the companion. 

We would like to thank to P. D. Sackett for making 
useful comments about the paper. This work was supported 
by the grant (20006-113-03-2) of the Korea Science & Engi- 
neering Foundation (KOSEF). 
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